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Abstract 

In this work we present a method to detect, identify and characterize stochastic 
information contained in an electrocardiogram (ECG). We assume, as it is weh 
known, that the ECG has information corresponding to many different processes 
related to the cardiac activity. We analyze scaling and Markov processes proper- 
ties of the detected stochastic information using the power spectrum of the ECG 
and the Fokker-Planck equation respectively. The detected stochastic information 
is then characterized by three measures. First, the slope of the power spectrum in a 
particular range of frequencies as a scaling parameter. Second, an empirical estima- 
tion of the drift and diffusion coefficients of the Fokker-Planck equation through the 
Kramers-Moyal coefficients which define the evolution of the probability distribution 
of the detected stochastic information. 
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1 Introduction 

The analysis of ECGs, or electrocardiography, remains the best non-invasive 
method for detecting or predicting coronary artery disease, it is also a very 
practical and nonexpensive method to obtain useful information about the 
quality of the cardiac activity. In the last few years, new concepts and methods 
from different fields of physics have been applied to the analysis of ECGs. 
The ECG can be considered as time series of measurements of an observable, 
i.e. the electrical activity of the heart recorded from the body surface. This 
approach has given some interesting results and provided new tools to define 
new measures of great potential for the characterization of the cardiac activity. 
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These new measures arc able to detect information contained in the ECGs 
not accessible with traditional methods. Such new useful information may 
be a mixture of both, stochastic and deterministic information, that has to 
be unraveled and analyzed with correspondingly different methods. In this 
work we present a method and some interesting results to analyze ECGs as 
time scries by detecting and characterizing the stochastic information in the 
context of Markov processes. In section 2 we present the problem of extracting 
stochastic information from a broad band power spectrum. In section 3 we 
present the method to characterize the detected stochastic information in the 
context of Markov processes applied to healthy and pathologic ECGs. Section 
5 is devoted to discussion and some conclusions. 



2 Detecting stochastic information in ECGs 

In a previous work we developed a method to separate deterministic and 
stochastic information from ECGs using the power spectrum, PS, of the ECG 
[1]. The fundamental idea was to improve the 1// behavior of the broad band 
power spectrum, BBPS, of the ECG in certain range of frequencies. This 
method uses an archetypal reconstruction of the ECG [2,3]. The difference 
between the original ECG and the archetypal reconstruction is considered the 
preliminary stochastic information, PSI. The final stochastic information, SI, 
was obtained by adjusting the power spectrum of the PSI to the best 1//" 
power spectrum. The 1//" power spectrum transformed back to the time 
space was considered the SI, (the difference between the 1//" and the power 
spectrum of the difference between the original ECG and its archetypal re- 
construction, transformed back to the time space, was considered part of the 
deterministic information). This process was partially successful because the 
detected SI obtained from the ECGs, showed interesting differences between 
healthy and pathologic ECGs when analyzed as Markov processes. However, 
the archetypal reconstruction is not necessarily a method to separate the de- 
terministic and stochastic information, it is a better or worse reconstruction 
of the whole ECG (depending on the size of the base in use) with its deter- 
ministic and stochastic components mixed. The reconstruction of the ECG 
with a scale dependent base, such as wavelets, complemented by a process of 
fitting the BBPS to a four parameter function is work in progress to improve 
the extraction of the SI from the ECG. In this work we use some preliminary 
results of the mentioned work in progress as the SI to be used and differentiate 
healthy and pathologic cases. Despite the preliminary character of the results 
used as the SI obtained from 21 healthy and 30 pathologic ECGs [4], we found 
important qualitative results of one measure that can distinguish different be- 
tween the healthy and pathologic cases where other measures do not capture 
such difference. This measure is the diffusion parameter of the Fokker-Planck 
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Fig. 1 Power spectrum of the ECGs. a) and b) correspond to the whole BBPS of six characteristic pathologic 
and six characteristic healthy cases respectively, c) corresponds to the same twelve BBPS in the frequency 
range ~ 7Hz to ~ 20Hz. 

equation estimated through the Kramers- Moyal coefficients and the method 
to estimate it is presented below with some results. 



3 Characterizing stochastic information from ECGs 

The ECGs have broad band power spectra, BBPS, with complex structure but 
with approximately 1//" behavior for an important range of the frequency do- 
main. The frequency range where the SI is going to be characterized goes from 
~ 7Hz to ~ 20Hz. Smaller and larger frequencies may be strongly affected by 
a variety of artifacts related with the different measurement and recording 
conditions for each ECG associated with nonstationarity and noise contami- 
nations respectively. Indeed, these two regions of the PS are more correlated 
to the data base used than to the healthy or pathologic character of the ECGs. 
The estimated value of a is a ffist characteristic of the SI obtained from the 
ECG. The different values of a estimated from all the SI obtained from each 
ECG do not permit us to distinguish between healthy and pathologic cases 
considering the mentioned range of frequencies. The following procedure will 
give us two more measures to characterize the SI obtained from from ECGs. 

These two measures are obtained by analyzing the stochastic information con- 
tained in the ECG in the context of Markov processes. First, we must give 
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Figure 2 Contour plots of the conditional probability distributions p(xi,Ti \ 2^3, T3) a) for the healthy case 
and b) for the pathologic case. 



evidences of Markov properties of the stochastic information. The following 
equation, Chapman-Kolmogorov, must hold if the time series under analysis 
corresponds to Markov processes [5]: 



where ri | X2, T2) is the conditional probability density of finding the value 
Xi at time ti given the value X2 at time T2. Evidence of such condition may be 
calculated from the time series and observed graphically in the contour lines 
of Fig. 2. Similarity of the left plot with the right one is evidence of Markov 
processes. 

The evidence of Markov process corresponding to the SI obtained from all the 
ECGs may be improved for certain values of the parameters of the method 
mentioned in section 2, compared with the archetypal reconstruction used 
in previous works to obtain the SI. With such evidences we can apply an 
empirical determination of the coefficients Di and D2 of the Fokker-Planck 
equation [7] 



p(Xi,ri I X3,T3) = / dX2p{x2,T2 \ X3, T3)p{xi, Ti \ X2, Tg) 
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Figure 3 Fokker-Planck coefficients, Di and D2 for the stochastic time series obtained from healthy a,b) 
and pathologic cases c,d). 

through the Kramers-Moyal coefficients 
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(3) 



where Dk{x,T) = limAT^oMk{x,T, At). 

The Fokker-Planck equation (2), defines the evolution of the probability dis- 
tributions p{x, t) of the Markov process at time scales defined by r. The 
coefficients Di and D2, drift and diffusion coefficients respectively, completely 
determine the Fokker-Planck equation [7,5,6]. Therefore, the coefficients Di 
and D2 can be used to characterize the stochastic information detected in the 
ECG as a Markov process. 

In figure 3 we present the estimates of the coefficients Di and D2 for the 
SI obtained from one representative pathologic ECG and one representative 
healthy ECG. We expect that the improvement of the procedure to obtain the 
SI and the quantification of these results will give a better and more roboust 
measure to distinguish healthy and pathologic ECG. However, in this work 
we report the potential of this measure and the method to obtain its values 
using the ECG as the unique source of information. 



4 Discussion and Conclusions 



Despite the preliminary character of the SI obtained from ECGs, the esti- 
mated values of the diffusion coefficient D2 allows as to identify interesting 
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new differences between healthy and pathologic cases. As it is observed in Fig. 
3, the plot of D2 as a function of x for the healthy cases did not present a 
global minimum for x = and a regular parabolic behavior as was the case for 
the pathologic cases. This observation is consistent with previous results and a 
general consensus that healthier cardiac activity corresponds to more complex 
and less regular behaviors, patterns, measures, etc. There were two exccep- 
tions to this observation, one pathologic and one healthy case. So far, we have 
not observed relevant information in the different estiamted values of the drift 
coefficient Di, that may be useful to distinguish healthy and pathologic ECGs. 
For all the ECGs analyzed, the corresponding BBPS can be qualitatively clas- 
sified in two distinctive groups corresponding more or less to the healthy and 
pathologic ECGs. However, the differences used to define these two groups 
correspond to very low and very high frequency ranges of the BBPS where 
external factors of the measuring and recording processes strongly affect the 
information contained in the EGG. In the frequency range between ~ 7Hz and 
~ 20Hz, we observe that all the BBPS are very similar qualitatively but very 
different quantitatively as indicated by the values of the measure a. However, 
the different values of a do not allow to make a clear clasification of healthy 
and pathologic cases. 
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